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We derive an approximate theory for Heisenberg spin ladders with two legs by mapping the spin 
dynamics onto the problem of hard-core 'bond-Bosons'. The parameters of the Bosonic Hamiltonian 
are obtained by matching anomalous Green's functions to Lanczos results and we find evidence for 
a strong renormalization due to quantum fluctuations. Various dynamical spin correlation functions 
are calculated and found to be in good agreement with Lanczos results. We then enlarge the effective 
Hamiltonian to describe the coupling of the Bosonic spin fluctuations to a single hole injected into 
the system and treat the hole-dynamics within the 'rainbow-diagram' approximation by Schmidt- 
Rink et. al. Theoretical predictions for the single hole spectral function are obtained and found to 
be in good agreement with Lanczos results. 
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I. INTRODUCTION 

Copper-oxides with a Cu02 plane containing line- 
defects, which results in ladder-like arrangmcnts of Cu- 
atoms 0, have recently received considerable attention 
1^. In addition to being interesting physical systems 
themselves [^,^, they may be considered, on the theoret- 
ical side, as an important stepping-stone to understand 
the fully 2D Cu02 planes of cuprate superconductors. 
Namely their special geometry makes two-legged ladders 
a realization of a particularly simple 'RVB'-type state, 
and developing a 'technology' for handling such RVB 
states clearly is one of the key issues in the theory of 
high-temperature superconductors. In the following, we 
want to derive a simple theoretical description of the spin 
dynamics and single-hole dynamics of ladders. Thereby 
we follow an approach which is similar in spirit to the 
Landau theory of Fermi liquids: we use a simple pic- 
ture of an 'RVB spin liquid' and its excitation spectrum 
and parameterize the dispersion of the spin excitations 
by few parameters, which are then obtained by matching 
the results to Lanczos calculations and finite-size anal- 
ysis. Using this parameterization we make quantitative 
predictions for various dynamical correlation functions 
which can be compared to exact diagonalization (and ex- 
periment). The theory actually allows to make rather de- 
tailed predictions, which in all cases studied are in good 
agreement with Lanczos results. We also discuss possi- 
bilities to extend the calculation to 2D systems. 
We consider the standard t — J model on a 2-leg ladder 
with N rungs. More precisely, the Hamiltonian reads 



H 



H.c. 



Si ■ Sj 



(1) 



Here denotes a summation over all pairs of near- 

est neighbors on the ladder, for bonds along the legs we 



choose ti 



t and Ji. 



t± and Jij = J±. 



- J, for bonds along the rungs 
In the following we adopt the 



values J± = J , t± = t, and J/t = 0.5, which may be 
roughly appropriate for actual the materials. We choose 
the X-axis along the direction of the legs, the y-axis along 
the rungs. 



II. SPIN DYNAMICS 

To begin with, we consider the case of half-filling and, 
following Ref. define the following operators: 



lj,X 

t 

t 



^(4T4a 4x4,?)' 
■^(4t4:T ~ 4j,4a) 

i 



:7^(4T4T + 4i4i)' 



(2) 



These create either a singlet or the three components of 
a triplet on the sites i and j. Following a large number 
of workers ||-|^, we will start out from the 'rung- RVB 
state', l^'o) — Yln=o 7i+y\^)' ^hich is the ground state 
in the limit J±/ J — > oo (here i + y denotes the nearest 
neighbor of i in y-direction) . Let us now consider the 
effect of switching on J. Obviously this will create 'fiuc- 
tuations' in the vacuum state, and it is of importance to 
clarify the nature of these. Due to the product nature 
of the vacuum it is sufficient to consider just one 4-site 
plaquette, see Figure ^ for the labelling of sites. For 



J=0 the ground state is |0) 



- „t J 



Sjslfac), with energy 



— (3/2)J^. It might appear that the fluctuation with 
the lowest cost in energy would be a '90-degree swap', 
i.e. a transition to a state with singlets along the legs: 
|l)=sj2'S43b^c). This state has energy — (3/2)J, i.e. it 
is degenerate with the vacuum in the limit J±=J. This 
line of thinking, however, is incorrect. The reason is that 
|1) is not orthogonal to the vacuum, more precisely one 
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finds (1|0) = 1/2. The most natural way to proceed is 
to form the orthogonal complement |2) = |1) — (1/2)|0), 
which (after normalization) is found to be 

The (manifestly rotation invariant) expression on the 
r.h.s. thereby is nothing but the Clebsch-Gordan combi- 
nation of two triplets along the rungs into a singlet. 




FIG. 1. Labelling of the sites in the 2x2 plaquette. 

This shows, that the 'true' excitation is the double 
triplet, rather than a state with singlets along the legs, 
or, put another way, that the state |1) is redundant and 
can be discarded if we keep |0) and |2). The same holds 
true for the last candidate state for a fluctuation, namely 
two triplets along the legs coupled to a singlet: 



|3) 



43, a: + ^12,j/^43,a 



+ ti2,zti3,z)\vac) 



Again we find that this state is not orthogonal to the vac- 
uum either, and orthogonalization again yields |2). The 
only fluctuation we have to take into account therefore is 
the formation of two triplets on adjacent rungs, the re- 
spective matrix element is {ti4.zt23,z\H\s\^sl^) — J/2. 
The energy increases by 2J±, which we interpret as 
twice the 'energy of formation' of a single triplet. Let 
us now assume that a double triplet has been formed, 
and consider its further development. First, the two 
triplets can 'recombine' and we return to the vacuum; 
the matrix element for this process again is J/2. Sec- 
ond, the two triplets can 'swap their species'. More pre- 
cisely, for aj^(3 there is a matrix element of the form 
(ii4,/3t23,/3|-ff|44,a43,a)=-'^/2- The last possibility is 
that one of the triplets interacts with a singlet on the 
as yet 'untouched' rung next to it. To study this in more 
detail we consider a 4-site plaquette containing spin 1 . It 
is convenient to form the states 



- —fit ct 

— 14,z''23 



± 4A3,z)\v<^c} 



which do have definite parities under the mirror opera- 
tions (Jx and (jy (see Figure ^). They are triplets with 
z-component 0. Next we consider which state could mix 



with one of these states. A first possibility would be two 
'rung triplets', coupled to a triplet: 



|4) = -^{tlA 



23, a 



ti4j23.x)\vac) 



However, this state has positive parity under x-reflection, 
whereas |±) both have negative parity. This state there- 
fore cannot be admixed [p|. The next possibility is 



t , 

12. z- 



43 



5l243,JI««c), 



(which does have the required negative a;-parity) , but this 
state is actually identical to |— ). Finally, one could think 
of two triplets along the legs coupled to form a triplet: 



|6) 



1 



T2' 



(*12,k43,v 42,v43,k)I^"'^) 



but, in fact, |6) = |+). The only process possible for 
the rung-triplet therefore is to 'propagate', i.e. exchange 
its position with a singlet on a neighboring rung. The 
hopping element is {{+\H\+) - {-\H-))/2 = J/2. The 
'90-degree swap' to a triplet along the leg (see the state 
1 5)) would give a redundant state and, most importantly, 
'anharmonic processes' whereby a triplet decays into two 
triplets (see the states |5) or |6)) are not possible either. 
Clearly, this absence of anharmonicity leads to a consid- 
erable simplification of the physics - it is the ultimate 
reason why, as will be seen later, the spin correlation 
function in a ladder is remarkably 'coherent'. 
There is one last process we need to discuss. It may hap- 
pen that two triplets of unlike species which have been 
created in different pair creation processes 'collide'. More 
precisely, for a ^ f3 there is a matrix element of the 
type (ti4_0i23.Q|J^Ki4 /s)"'^/^' which describes two 
triplets of unlike species 'hopping over' one another. 
We now represent the presence of a triplet on the n*'' 
rung by the presence of a 'book-keeping boson' created 
by t}^ ^ . Grouping the three possible triplet states into 
a single 3- vector, tn—{tn,x,tn,y,tn,z), we can describe all 
processes involving non-redundant states by the following 
manifestly rotation-invariant Hamiltonian 

H = J^Y.*n-tn 

n 

+ ^ E( *" • + H.c. ) + ^ 5]( 4 • K+i + H.c. ) 

n n 
n 

This Hamiltonian has also been derived by Gopalan et al. 
Q . The form of the Hamiltonian (^ suggests a quite ob- 
vious approximation, namely to break down the quartic 
terms in a BCS-like fashion, e.g.: 



^ — (^n+l.Q^n,a)- 



(4) 
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As a next step, we might assume that the hopping inte- 
grals and pair creation matrix elements in (j^) are renor- 
malized in a Gutzwiller-like fashion, to mimic the effect 
of the hard-core constraint. After Fourier transform, we 
would thus arrive at an approximate Hamiltonian of the 
form 



fc>0 



-iAkti-tl, + H.c.)]. (5) 



Assuming that the t in this Hamiltonian are free Bosons 
this is readily solved by the ansatz 

7Li,,„ ^Wfetfe + Ufetlfe, (6) 
to give the (3-fold degenerate) dispersion 



(7) 



The parameters Ck, should be calculated in some 
approximate fashion, by mean-field and Gutzwiller- 
approximation to the original Hamiltonian j5|. We 
note, however, that there are a number of difficulties with 
such an approach: the density of Bosons nf, can be ob- 
tained from the spin correlation function along the rungs 
using the identity —3(1 — ri;,)/4-|-ri(,/4 = (Si-Si+y). The 
spin correlation function thereby can be estimated from 
exact diagonalization results. For Jj_/J= 1 we thus ob- 
tain rib ~ 0.3. The quantum fluctuations thus are strong, 
and we may expect that the matrix elements for propa- 
gation and pair creation in (||) are heavily renormalized. 
On the other hand, the matrix elements of the quartic 
terms in (^) will not be renormalized at all due to the 
hard-core constraint, because these terms do not change 
the Boson occupation of any site. We may thus expect a 
subtle interplay of Gutzwiller projection and mean-field 
decomposition and to avoid poorly controlled approxi- 
mations we resort to numerical techniques to extract the 
parameters of the Hamiltonian (^ and, in doing so, more- 
over check the quality of the mapping to noninteracting 
Bosons per se. To that end we define the operator 

If the rung (n, n + y) is in a singlet state, transforms it 
into the 2-component of the triplet; if the rung is in any of 
the triplet states, the state is annihilated. This operator, 
while acting entirely within the Hilbert space of the orig- 
inal ladder system, thus may be viewed as a realization 
of the hard-core Boson creation operator. Next, using 
standard Lanczos techniques, we evaluate the following 
Green's functions: 



AM(k,uj) ^^-('^o\tI 7 ^-^ ^T^I*o), (8) 
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where |5'o) (-Eo) denote the ground state wave function 
(energy) of the ladder and is the (1 dimensional) 
Fourier transform of t}^. Our goal is to map the spin 
excitations of the ladder onto a system of free 'Quasi- 
Bosons' goverened by the Hamiltonian (||). If we want to 
compare spectral properties, we have to take into account 
that the spectral weight of the hard-core Bosons may be 
strongly renormalized. For example, as a rigorous iden- 
tity we have ([r.;, t/]) = 1 — 4nf,/3R:!0.6, rather than 1 as it 
would be for free Bosons. To discuss the spectral func- 
tions we therefore assume that tj,— s-'v/Ztfe, where t^ is a 
free Boson operator, and Z the wave function renormal- 
ization constant. 

Assuming then that the free Bosons tk are indeed de- 
scribed by a Hamiltonian of the form (|^) one can derive 
the following expressions for the above Green's functions: 



w(A_|_(fc, uj) 



L0Air,t{k,Uj) 

- A_(fc, w) 



ZekS{uj - ujk) 
Z6{uj - LOk). 



(9) 



In other words, the wave function renormalization Z as 
well as the energy and pair creation amplitude can 
be read off from the dispersion of the pole strength in the 
Green's functions (||). Then, using the obtained values 
of tk and Afc to calculate the dispersion of the excitation 
energy from (|^) and comparing with the actual numeri- 
cal values should provide a stringent cross-check for the 
validity of the mapping to the free-Boson Hamiltonian 
(^. Figure ^ compares Aj^{k,uj) and uj A^ntik , lij) and 
uj{A+{k,uj) -f A_(fc,w)). To begin with, the A+{k,uj) 
shows a series of sharp dispersive peaks, which how- 
ever trace out a rather unusual dispersion with a shal- 
low local minimum at fc=0. Next, lo Aint{k ^ to) has a 
pronounced low energy peak too, which coincides with 
that of A^{k,oj). The residuum changes sign at approxi- 
mately fc = 7r/2 (note that Aint is an off-diagonal matrix 
element of the resolvent operator; its residuum therefore 
need not be positive definite), the magnitude is nearly 
identical at k~Q and k—i:. This suggests a fc-dependence 
of the form A^, cx cos(fc), as one would expect for a 
pair amplitude corresponding to pair creation on near- 
est neighbors. Similarly, Lj{A^{k,uj) + A^{k,uj)) shows 
low energy peaks coniciding with those of A_|_(fc, uj). The 
dispersion of the weight, however, is more complicated 
than for the pairing amplitude. 



3 




FIG. 2. Spectral functions A+{k,Lu) (full line), 
Lj{A+{k,u)) + A^{k,u>)) (dashed line) and u)Aint{k,u}) (dotted 
line) obtained by Lanczos diagonalization of a 2 x 10 ladder. 
(5- functions are replaced by Lorentzians of width O.IJ. 



For a quantitative analysis, we proceed to Figure |^. 
Part (a) shows first of all the weight Z obtained from 
A+{k,uj) — A^{k,uj). Remarkably enough, Z is quite in- 
dependent of k and moreover close to the estimate of 0.6 
obtained from ([t^jt/]). The fact that Z is indeed fairly 
fc-independent is a first indication for the applicability 
of a free Boson Hamiltonian. To avoid a proliferation 
of adjustable parameters we take the fc-average, which 
is Z=0.53 (nearly independent of N for > 6). Using 
this value of Z we then calculated the ik and and 
Fourier transform them with respect to k. For simplic- 
ity, we terminate the Fourier series after the third (i.e. 
cos(2fc)-like) term. This gives already a very satisfactory 
fit to the data, as seen in Figure |[ Figure ||b shows 
the comparison of the numerical excitation energy and 
the dispersion calculated from (^, using the values of 
ik and Afe extracted from the weights of the correlation 
functions. There is quite reasonable agreement, which 
indicates that the effective free-Boson Hamiltonian itself 
is a quite good description of the spin dynamics. 
We proceed to extrapolate the results to the infinite 
chain. To that end, we performed the calculation for 
different N and calculated the Fourier coefficients of ik 
and Afc. Then, Figures ^(c) and (d) show plots of these 
Fourier coefhcients vs. 1/N'^. The plots rather obviously 
suggest that all of the Fourier coefficients can be 



2.5 



1.5 - 



0.5 



-0.5 



0.4 
0.3 
0.2 
0.1 



■0.1 
-0.2 
-0.3 

-0.4 



.111111. 


.111111. 






-1 1 1 1 1 ""l- 


- 

r ~ 

-I 1 1 1 1 i- 


0.2 0.4 0.6 0.8 1 

k/7T 


0.2 0.4 0.6 0.8 1 

k/7i 


1 ' 1 ' 1 ' 

(c)- 


- ■ ' ' (d)- 








- 


— ■ 


€„-2 " 












■ ^ 


" , 1 , 1 , 1 , " 


" , 1 , 1 , 1 , " 



0.01 



0.03 



0.02 0.03 0.01 0.02 

FIG. 3. (a): (triangles), (squares), and Z (hexagons) 
as obtained from the pole-strengths of the Greens functions 
plotted versus momentum. The continuous line gives Fourier 
expansions with the lowest 3 harmonics. 

(b) : ujk calculated from (Q) (triangles) compared to the exact 
excitation energies (squares). 

(c) and (d): Scaling of the Fourier coefficients of A and e with 
the length of the ladder. 



expressed as c{N) = Cqo + c' /N"^ to good accuracy. The 
extrapolated values are given in Table |. We note in 
passing that these values give a 'spin gap' (i.e. the energy 
of the triplet with fc = 7r)of0.51J, in agreement with ex- 
act diagonalization |9| and DMRG calculations Let 
us now briefly discuss these values. For the pairing ampli- 
tude Afc we have a by far dominant cos(fc)-component, 
Ai, consistent with pair creation/annihilation on near- 
est neighbors. The second- nearest neighbor amplitude is 
substantially smaller, and also the uniform component is 
very small. 

The results for the rcnormalized energy are more sur- 
prising. To begin with, the constant term is 1.77 J , rather 
than J±—J, as one would expect. A possible explanation 
for this strong increase is that a triplet on (say) the rung 
number i blocks the 'pair creation' of triplets on the pairs 
(i — and + 1). The corresponding loss of fluc- 
tuation energy increases the energy cost for creating a 
triplet, i.e. the on-site energy of a triplet. In addition, 
presence of a triplet will partially inhibit the motion of 
other triplets and thus cause a further loss of delocaliza- 
tion energy. Quite obviously these effects combined are 
quite substantial. Next, the coefficient of the nearest- 
neighbor harmonic is rather small (only 0.21). This may 
be understood simply in terms of a Gutzwillcr-like re- 
duction of the mobility due to the 'excluded volume' 
occupied by other triplets. Moreover, the terms which 
describe the 'exchange hopping' and effectively propa- 
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gate a triplet by one lattice site may interfere destruc- 
tively with the 'ordinary' motion of a triplet, and thus re- 
duce the effective hopping integral. A further surprise is 
the large amplitude for second-nearest neighbor hopping. 
Here one could envisage processes where the propagating 
triplet encounters a pair of triplets, recombines with one 
of them and thereby transfers its momentum to the re- 
maining triplet, so that the propagating triplet effectively 
has been transferred by two sites. Such a process would 
be proportional to the probability of finding a quantum 
fluctuation, which is rather high in the present case. All 
in all, the data show that the propagation of the triplets is 
strongly renormalized in the relatively dense and strongly 
interacting 'background' of the other triplets. 
Assuming that the 'effective Hamiltonian' (||) with the 
extrapolated parameters in Table | gives a good descrip- 
tion of the spin dynamics, we proceed to calculate dy- 
namical correlation functions relevant to experiment. We 
consider the spin operator on a single rung and 'trans- 
late it' into the language of the hard-core Bosons. The 
operator 5'Q{A;_L=7r) turns a singlet into an a-triplet and 
vice versa; the other two types of triplets are annihi- 
lated, whence: S{k±=TT) — > t''+t. Next, the operator 
'S'z(fc_L=0) annihilates a singlet and converts e.g. t\.~^it'^y 
(compare (||)). S{k±=Q) therefore must be bilinear in 
the triplet operators and the only possibility is t'^ x t. 
From its acting on t|, the prefactor must be —i, whence: 
S{k^—Q)^—it'^ X t. Obviously this is also the only way 
to construct a Hermitean operator (5(fcx=0) is the op- 
erator of total spin on one rung) from the triplets. A 
subtle point is the renormalization of the spectral weight: 
whereas S^(A;_L=7r) actually changes the number of Quasi- 
Bosons and thus should be renormalized in a similar fash- 
ion as the hard-core Boson addition and removal spectra, 
S(fc_L=0) does not change the Boson occupation of any 
rung and thus should remain essentially unrenormalized. 
Based on these considerations, we multiply S{k±=Ti) by 
but leave S{k^—Q) unrenormalized. Then, upon 
Fourier transformation of the single rung operators and 
switching to the bond-Bosons we obtain the spin corre- 
lation functions: 



S{q,-K,z) = 



Z- LOk 

[Uk+qVk - U^kV-{k+q)Y 
Z - {uJk+q + UJk) 



(10) 



Note that S{q, 0, z) ^ for g ^ 0; this is what must 
come out because unlike conventional spin-wave theory 
the 'rung-RVB' state is an exact singlet, and both Hamil- 
tonians, (^ and (j^) are rotationally invariant; the ground 
state thus is an exact singlet, whence the operator of to- 
tal spin must annihilate it. The obtained spin correlation 
functions are shown in Figure ^. We first note the rather 
different intensity of the two correlations functions. As 
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FIG. 4. Spin correl&tion function ^(q, a/J- i ■ O.IJ) obtained 
by numerical evaluation of (|l^ (full line) compared to the re- 
sult of Lanczos diagonalization on an N=12 ladder (dotted 
line). Spectra for qy—O are multiplied by a factor of 5. 

could be expected on the basis of ([To|), the correlation 
function for qy=iT consists essentially of a single peak, 
whereas the spectrum for qy=0 has more 'cusp-like' ap- 
pearance. It is interesting to note that these differences 
are reproduced quite well by the results of Lanczos cal- 
culation, which are also shown in Figure (Q): the qy~TT 
spectrum is indeed remarkably sharp, with nearly all the 
spectral weight concentrated in just one peak for every 
momentum. It should be noted that the agreement of the 
dispersion is not really surprising, this was actually seen 
already in Figure ^. The dispersion of the peak intensity, 
however, also agrees very well with the Lanczos result 
and thus provides further evidence for the correctness of 
the mapping to the Boson Hamiltonian. 
By contrast, the spectra for qy = usually consist of 
several peaks and one can already envisage how these 
spectra develop into the cusp-like spectra produced by 
our theory. The dispersion of the 'peaks' for qy—0 is in 
good agreement with theory, showing a shallow maxi- 
mum at qx~3Tr/5 and approaching zero for q^ — > 0. All 
in all, the spin correlation function is obviously very well 
reproduced by our calculation. We proceed to a discus- 
sion of the 'energy correlation function' ea{q,oj), i.e. the 
dynamical correlation function of the operator 



Si • s 



i+6i • 



Apart from being a further probe of the spin dynamics, 
this correlation function plays a key role in the theory 
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of phonon-assisted two-magnon absorption by Lorenzana 
and Sawatzky [0. It therefore can be probed experi- 
mentally in infra-red absorption measurements. Here we 
restrict ourselves to the operator Hy{q), which is easy to 
'translate' to the Quasi-Boson system. Namely one can 
replace 

Si ■ Si+a > t| • t j 

whence 

, , _ 3 {uk+qVk + U-kV-(k+q))'^ 

Since Hy{q) does not change the number of Bosons on 
any rung, we have not added a factor of Z . The other 
correlation function, ex{q,oj) is more difficult. In princi- 
ple, the exchange along a given bond along the legs of the 
ladder is nothing but the respective part of the effective 
Hamiltonian (|^). However, as has been discussed above, 
if we go over to the simplified Hamiltonian (|^) there is 
quite a strong renormalization and also a generation of 
additional next-nearest-neighbor hopping terms. We be- 
lieve this is difficult to treat in a reasonably controlled 
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FIG. 5. Energy correlation "Wnction e{q,u) — i ■ O.IJ) ob- 
tained by numerical evaluation of (dotted line) compared 
to the result of Lanczos diagonalization on a A/^ = 12 ladder 
(full line). Note that Qy — n cannot be meaningfully defined 
and is therefore omitted. 

way and therefore omit a discussion of this correlation 
function. Then, ey{q,Lu) is shown in Figure || and com- 
pared to Lanczos results. Unlike the spin correlation 
function, the detailed agreement of the lineshape between 



the theory and Lanczos is not very good, with the dis- 
crepancy being particularly large around q=0. There, 
the Lanczos result shows an intense peak at an energy 
of w 2.5J, which is completely absent in the theoreti- 
cal spectra. We believe that this peak corresponds to a 
'bi-triplet', i.e. a bound state of two triplets on nearest 
neighbors. Two triplets on nearest neighbors can effi- 
ciently lower their energy by using the exchange-like in- 
teraction described by the quartic terms in (^). This 
may produce a virtual bound state similar as the bi- 
magnon excitation familiar from the Raman spectrum of 
the 2D Heisenberg antiferromagnet. Our simple particle- 
hole picture for the correlation function naturally does 
not take into account such a bound state formation and 
therefore cannot reproduce this resonance. On the other 
hand, for momenta which are more distant from q—0 the 
agreement between Lanczos and theory improves some- 
what. The theoretical spectra reproduce the lower edge 
and 'width' of the spetra reasonably well, and also the q- 
dependence of the total spectral weight is approximately 
correct. 



III. HOLE DYNAMICS AND PHOTOEMISSION 

Next, we proceed to a theory for photoemission [^ . 
To that end, we need to study the motion of a single hole 
in the 'spin background'. The previous discussion has 
shown that we needed to retain only singlets and triplets 
along the rungs and that, as far as the spin dynamics is 
concerned, all other possible states are redundant. This 
suggests to discuss the added hole in such a way that 
it 'fits' with the rung basis used for discussing the spin 
dynamics. We define bonding and antibonding states of 
one electron along a rung: 

^lk,,a = ^(^l<^ + ^'''^i+y,.)' (11) 

where we have represented the parity under as mo- 
mentum or TT in y-direction. These states have an 'on 
site energy' of — cos{ky)t± and we now consider their 
propagation. The goal thereby is to interpret a singly oc- 
cupied rung as being occupied by an 'effective Fermion', 
which has the z-spin and ky of the respective single- 
electron state. Then we want to enlarge the effective 
Hamiltonian ^ by terms describing the propagation of 
these effective Fermions as well as their interaction with 
the triplet excitations. 

The form of the terms in this effective Hamiltonian can 
be inferred already from the requirements of positive par- 
ity under a^, spin rotation and time-reversal invariance 
and Hermiticity. We note that the rung-singlet and a|j q 
have positive parity under a^, whereas the triplets and 
aliT^ „ have negative parity. Then, from the spinors a we 
can construct two 3-vectors 
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n.ky-\-7T,(T' : 



(12) 



where r is the vector of PauU matrices. Thereby Sn.m 
has odd parity under exchange of the legs, Sn,m has even 
parity, and we have Sl^ ^=Sm,n- Both vectors are odd 
under time reversal (as is t). The processes we need to 
describe are shown schematically in Figure O: by virtue 





C3 ( 





C4 



FIG. 6. Possible interactions between a hole and a doubly 
occupied rung. Each process is labeled by the prefactor of 
the terms in the Hamiltonian (|l^) which describes it. The 
possible processes are: exchange with a singlet (ci), exchange 
with a singlet which is transformed into a triplet (C2), ex- 
change with a triplet whereby the triplet remains unchanged 
(ca), exchange with a triplet whereby the triplet changes its 

(C4). 

the hopping term a singly occupied rung can exchange 
its position with a 'fully occupied' rung. Thereby the 
fully occupied rung can remain in a singlet, change from 
singlet to triplet (and vice versa), remain in a given 
triplet state, or remain in a triplet state but change its 
Sz- These processes can be described by the following 
coupling terms, whose form follows simply from the re- 
quirements of spin rotation invariance, even parity, and 
Hermiticity: 



C2 ^[ til • {riji+lSn.n+l + nn-lSn,n-l) + H.C.] 



+ C3 • *" H al,ky,aan-Uky,a + H.C. ] 



n ky ,{T 

C4^[iS'„,„-i • (tj'j X t„_i) +H.C. ], 



(13) 



where 



Tlji — 1 • tji . 



To determine the numerical values of the coefficients 
we start with the state |1) = a\ f. |sj|uac). Denoting the 
hopping term along the legs by fit we find: 



Ht\l) = 2(4,fc„,T4+4,fc„+7r,T4,2-v^4,fc„+,ra4,+)bac). 

Next, starting from the state |1') = a| -i^tX z\''^a,c) we 
obtain 

Ht\l') = ^(a2,fc„+^,T*^+4,fc„T4,^ + ^/2a^^^^4,+)kac). 

Using these two equations we find ci=</2, C2^t, C3=</2, 
and C4=— t. 

Next, we consider the action of the exchange terms along 
the legs. Namely the electron on a singly occupied rung 
can exchange with an electron belonging to a singlet or 
triplet on a neighboring rung. Again we can write down 
the general form of the Hamiltonian as 

H^c',J2[ {ti+, + tLi) • 5„,„ + H.c] 

n 
n 

Denoting the exchange term along the legs as Hj, we find 

^j|l> = ^ialky+^,A,z - ^«lfc„+7r,i4,+)l«ac), 

Hj\l') = ^{aiky+.A4 + V2a^^,^4,+)l««c), (15) 

whence C2= J/2 and 04=— J/2. Together with the Hamil- 
tonian for the triplet-dynamics, Eq. (^) then gives a com- 
plete description for the hole motion in the ladder. 
Having found all possible ways of interaction between the 
hole and the spin excitations, we need to consider a suit- 
able approximation to handle these terms. To begin with, 
by treating the terms oc ci , C2 in ( p^ ) in a mean-field like 
way, we may hope to obtain a 'renormalized hopping' in- 
tegral for the hole. Namely we replace n„ — > (1 — n-i,), 
with rib = 0.3 the density of triplets (see above). Next, 
in the term oc C3 we replace 



''n-l 



-^^cos(fc) vl. 



(16) 



We thus obtain the 'effective hopping integral' 

k 

Numerical evaluation shows that this is a very small 
quantity, ie//=0.27i, the reason being that the second 
term on the r.h.s. is negative and of quite appreciable 
magnitude. 

The form of the terms cx C2, C2 suggests the 'rainbow di- 
agram' approximation due to Schmidt-Rink et al. |l3[ to 
treat them. This still leaves us with the terms (x 04,04; 
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detailed investigation shows, however, that the matrix 
elements of these terms are proportional to higher pow- 
ers of the (small) coherence factor Vk; we take this as a 
justification to neglect these terms. 

Next, we briefly discuss additional 'renormalizations' 
which might occur. In the same way as a triplet on some 
given rung will 'block' quantum fluctuations and hamper 
the propagation of other triplets, a singly occupied rung 
will do the same. For the triplets, this effect has led to 
a quite dramatic increase of the 'energy of formation' of 
the triplet. However, the corresponding renormalization 
of the on-site energy for the hole-like Fermions is most 
likely independent of the ky of the hole, and hence can 
be absorbed into an overall constant shift of the spectral 
function. Since, the hole will not be able to propagate 
by two lattice sites by coupling to a quantum fluctua- 
tion, so that we expect that unlike the case of spin ex- 
citations there will be no 'dynamically generated' next- 
nearest neighbor hopping. 

Performing the Fourier and Bogoliubov transform wc fi- 
nally arrive at the following Hamiltonian 

k.ij q 

efc 2te// cos(fc^) - <_L cos(fcy), 
m{k,q) = 5q^,^ [ {2tcos{kx - Qx) + J cos{qx))uq^ 

-{2tcos{kx) + Jcos{qx))vq^ ]. (17) 

Thereby it is understood that ujq has only one branch 
with Qy—O. Formally, this is already very similar to the 
Hamiltonian derived by Schmidt- Rink et al. [l3| for hole 
motion in a Heisenberg antiferromagnet; the only dif- 
ferences are that we have spinful holes, a nonvanishing 
dispersion Ck already for the 'bare holes', and 3 branches 
of spin excitations rather than 1. In short, the Hamilto- 
nian is explicitely spin-rotation invariant, as it has to be 
in a 'spin liquid'. Technically this does not make any dif- 
ference and the equation for the self-energy I](fc, lu) reads 
now 



q 



\m{k,q) 



i;(fe -q,uj)' 



(18) 



The self-consistency equation can be solved for relatively 
large systems (we used N =200), and the self-energy be 
used to calculate the hole-like photoemission spectrum 



A{k,uj) 



TT - 



1 



efc-S(fc,w)-)-iO+' 



(19) 



If we want to compare this spectrum with Lanczos re- 
sults some care is necessary. The reason is that 
the spectral function we have calculated above is the one 
for the creation of a 'bare' hole. When calculating the 



photoemission spectrum, i.e. the spectral function of the 
operator Ck^m there exists the possibility that the anni- 
hilation operator 'hits' a rung in a triplet state. This will 
lead to terms of the form at , „ 



''k+q<j''q photoemission 
operator, i.e. the creation of the hole is accompanied by 
the annihilation of a spin excitation. If we want to have 
numerical results for the bare hole spectral function, we 
therefore should use the Fourier transform of the opera- 
tor 



+e'^y{ c„+g,tn„4 - Cn+y^iS~ ) 



(20) 



This operator replaces a rung-singlet by a singly occu- 
pied bond with the proper ky, and annihilates any triplet. 
This may therefore be considered as a creation operator 
for a 'bare' hole. Then, Figure (^ compares the 'bare 
hole' spectral function obtained by the self-consistent 
Born approximation and the Lanczos spectra of the op- 
erator Ck^a- The agreement is obviously quite good, al- 
though the self-consistent Born result tends to produce 
'too coherent' spectra and does not put enough weight 



1 1 1 1 1 1 1 1 1 1 1 1 

k=(Tr,0) 




' 1 ' 1 ' 1 ' 
k = (Tr,7r) 


1 ' 1 ' 1 ' 


k=(47r/5,0) 




i 




k = (37r/5,0) 




k = (37l/5,7T) 




k = (2n/5,0) 


k = (87l/5,7T) 


k = (n/5.0) 


k = (7r/5,7i) 


k = (0,0) 

1 1 .l-^l\-_l-^-l 'l<^f^ , 


k = (0,7T) 





-5 -4 -3 -2-10 1 -5 -4 -3 -2-10 1 
FIG. 7. Hole sp^ttral function by se4f-consistent Born 
aproximation (full line) compared to the spectrum of Cfe_o- ob- 
tained by Lanczos diagonalization in a 2 x 10 ladder (dotted 
line). The ratio J/t=0.5. 

into incoherent continua. This may be an indication 
that our renormalization of the coupling matrix element 
TO(fc, q) by (1 — rib) is too strong - a larger coupling 
would presumably lead to more incoherent weight. The 
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key features, however, are essentially identical in both 
spectra: the intense band in the ky=0 sector, which dis- 
perses slightly upwards from k=0 towards k—TT/2 and 
then more or less levels off; the equally intense band in 
the ky = TT sector, which starts out at fc=7r, disperses to- 
wards lower energy and quickly becomes 'overdamped'. 
The SCB results somewhat underestimate the bandwidth 
of the ky—0 quasiparticle band and put the ky=0 at a 
somewhat too high energy. The former could probably 
be remedied by adjustment of the 'bare hole hopping in- 
tegral' teff - however, since there is no rigorous way to do 
so, we decided not to do this. Apart from these relatively 
minor problems, however, there is quite good agreement. 
As a last remark we note that a comparison with the 
'full' photoemission spectrum jl2| shows quite substan- 
tial differences in the spectral weight of some features. 
The additional processes where the photoemission oper- 
ator annihilates a quantum fluctuation (in this case a 
triplet) thus are quite important for a discussion of the 
true photoemission lineshape (see also Refs. for a 

discussion of this issue in the 2D systems). 

IV. CONCLUSION 

In summary, we have derived a simple theory of two- 
legged spin ladders which reproduces a number of nu- 
merical results quite well. Being of comparable simplic- 
ity as linear spin wave theory for the planar Heisenberg 
antiferromagnet, the theory nevertheless allows to make 
quantitative calculations of physical quantities, which in 
all cases compare favourably with the results of Lanczos 
diagonalization. Just as linear spin wave theory may be 
viewed as constructing an effective Hamiltonian for the 
pair creation and propagation of fluctuations around the 
'Neel- vacuum', the present theory may be viewed as an 
entirely analogous expansion around the 'RVB-vacuum': 
the Hamiltonians (||) and (^) may be thought of as de- 
scribing the dynamics of triplet-like fluctuations around 
the RVB vacuum, and to discuss the hole dynamics we 
basically needed to describe the coupling between these 
triplet-excitations and the doped hole. The great sim- 
plification which made the calculation possible was the 
special geometry of the 2-legged ladder which immedi- 
ately suggested a unique RVB-vacuum around which we 
could 'expand' the fluctuating ground state. We note 
that there actually exists another system with such a 
unique RVB-vacuum, namely the Kondo lattice. For this 
system an analogous procedure involving Fermionic fluc- 
tuations is possible and leads to excellent results when 
compared to Lanczos diagonalization [l^ . 
In a two dimensional t — J-like system, such a unique 
vacuum docs in general not exist, unless one were to as- 
sume some spatial inhomogeneities such as stripes [ pT[ . 
Rather, for a translationally invariant state, the most 
natural vacuum would be a statistical average over all 



possible dimer-coverings of the plane, where each dimer 
corresponds to a singlet . The analogue of the triplet- 
like excitations in the ladder then would be states where 
one or several singlet-dimers are substituted by a triplet, 
and this excited dimer propagates through the system. In 
fact, up to an additional statistical averaging to account 
for the probability of suitable dimer configurations, the 
calculation of matrix elements for the pair creation and 
propagation of these excited dimers proceeds in an en- 
tirely analogous way as for the ladder . We note that 
such a picture for the low energy excitations of the RVB 
state is almost mandatory in the framework of the SO {5) 
symmetric theory of cuprate superconductors by Zhang 
: there, the 7r-operator, which actually accomplishes 
the SO (5) rotations of spin-excited, half-filled states into 
doped, superconducting ground states precisely con- 
verts excited dimers with momentum (7r,7r) into d^2_y2 
symmetric hole pairs with momentum (0,0). Then, the 
simplest picture of the SO (5) rotation from the antiferro- 
magnetic to the superconducting state would be that 
a condensate of excited dimers with momentum (tt, tt) at 
half-filling (i.e. a state with nonvanishing staggered mag- 
netization) is converted into a condensate of c?-like hole 
pairs (i.e. a c?-wave superconductor) in the doped case. 
A description of the spin excitations in 2D in terms of 
Boson-like excited dimers as for the t — J ladder thus 
may be a very natural starting point for this promising 
scenario. 

Instructive discussions with Dr O. P. Sushkov and Dr E. 
Arrigoni are most gratefully acknowledged. 



TABLE I. The Fourier coefficients of the renormalized en- 
ergy Efc and pairing amplitude A^, extrapolated to infinite 
length of the ladder. 
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1.77 


0.21 


-0.32 


Afe 


0.07 


1.03 


-0.17 
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